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ABSTRACT 



Context. When modeling infrared or y-ray data as a linear combination of observed gas tracers, excess emission has been detected 
compared to expectations from known neutral and molecular gas traced by HI and CO measurements, respectively. This excess might 
correspond to additional gas component. This so-called "dark gas" (DC) has been observed in our Galaxy, as well as the Magellanic 
Clouds. 

Aims. For the first time, we investigate the correlation between visible extinction (Ay) data and gas tracers on large scales in the solar 
neighborhood, to detect DG and to verify our compatibility with previous studies. 

Methods. Our work focuses on both the solar neighborhood (|b| > 10°), and the inner and outer Galaxy, as well as on four individual 
regions: Taurus, Orion, Cepheus-Polaris, and Aquila-Ophiuchus. Thanks to the recent production of an all-sky Ay map, we first 
perform the correlation between A v and both HI and CO emission over the most diffuse regions (with low-to-intermediate gas 
column densities), to derive the optimal (Ay/NH) lcf ratio. We then iterate the analysis over the entire regions (including low and high 
gas column densities) to estimate the CO-to-H 2 conversion factor, as well as the DG mass fraction. 

Results. The average extinction to gas column-density ratio in the solar neig hborhood is found to be (A v /N H ) ref = 6.53 x 10~ 22 
mag. cm 2 , with significant differences between the inner and outer Galaxy, of about 60%. We derive an average value of the CO- 
to-H2 conversion factor of Xco = 1-67 x 10 20 H2 cirT 2 /(K km s -1 ), with significant variations between nearby clouds. In the solar 
neighborhood, the gas mass in the dark component is found to be 19% relative to that in the atomic component and 164% relative to 
the one traced by CO. These results are compatible with a recent analysis of Planck data within the uncertainties of our measurements. 
We estimate the ratio of dark gas to total molecular gas to be 0.62 in the solar neighborhood. The HI-to-H 2 and H 2 -to-CO transitions 
appear for A v - 0.2 mag and Ay — 1.5 mag, respectively, in agreement with theoretical models of dark-H 2 gas. 

Key words. ISM:dust, extinction - ISM: clouds - solar neighborhood 



1. Introduction 

The interstellar medium (ISM) is composed of gas and dust, 
where the dust component represents about 1% of the total 
ISM mass. Interstellar gas provides the material for star for- 
mation, making it a central part in the life cycle of matter. 
The gas is either atomic, molecular, or ionized and these 
three phases are commonly observed using the HI 21 cm 
emission, CO transitions, and either H a or free-free emission, 
respectively. Several studies have reported an emission excess 
in the far-infrared (FIR) with respect to the gas, that is cor- 
related with none of the gas tracers. Thes e excess det e ctions 
were first obtained for G a lactic regions (Rea ch et al 
Meverdie rks & Heithausenl 11996). By analyzing the diffuse 
y-ray emission from the ISM, iGrenier et al.l (12005b concluded 
that there is an additional gas phase in our Galaxy called "dark 
gas" (DG), which has a non-negligible mass. Thi s component 
has also been discovered in the Magellanic clouds dLerov et all 
I2007L iBernard et all 120081: iRoman -Duval et al .1 120101) . 



The iPlanck co llaboration (2 01 lal) obtained all-sky maps 
of the dust optical depth and, comparing it with the observed 
gas column density, constructed a map of the DG distribution 
covering a large fraction of the intermediate Galactic latitude 
sky (|b| > 10°). On average, they estimated that the mass of 



DG is 28% of the atomic mass and 118% of the molecular 
gas traced by CO emission in the solar neighborhood. These 
results indicate that the DG is detected at intermediate hy- 
drogen column densities, corresponding to extinctions values 
between 0.4 and 2.5 mag. The first possible explanation is 
that the DG could be molecular gas, which is not detected 
using the CO(J = 1-0) transition. It is indeed expected that a 
layer of pure H2 or dissoc i ated C O should exist around dense 
clouds (e.g. IWolfire et all [2.010) that would be undetectable 
because at temperatures below 100 K (which is typical of these 
environments), the fluxes from the H2 rotational transitions are 
too low. However, this la yer could als o be tr aced using other 
species such as C + (see Lange r~et all l2010[). Other possible 
origins were suggested by the Planck co llaboration! (1201 lab 
to explain the observed departure from linearity between the 
dust optical depth (t) and the observable gas column density, 
namely: (1) variations in the dust/gas ratio (D/G), (2) weak 
CO emission undetected at the sensitivity of the CO survey, 
(3) optically thin approximation for the HI emission not valid 
over the entire sky, and (4) formation of dust aggregates, in- 
ducing a higher dust emissivity. For hypothesis (1), the authors 
concluded that variations in D/G of about 30% are unlikely in 
the solar neighborhood. Assumptions (2) and (3) were tested 
by performing the analysis with an upper limit in the weak CO 
emission in one case, and with a different HI spin temperature 
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in the other case. The authors deduced that these two effects 
could not account for the whole excess. The last option (4) 
can be tested using extinction data, since dust aggregates are 
expected to have a higher FIR emissivity than isolated grains, 
but mostly unaffected optical properties in the visible and 
ultra-violet. Therefore, no substantial Ay excess relative to the 
gas column density should be observed in that case. We note 
that y-ray observations would not detect any excess if the FIR 
excess were due to either D/G or dust optical property variations. 

The large dust particles that emit in the FIR are also re- 
sponsible for the main extinction in the visible and near-infrared 
(NIR). Extinction data are therefore well-suited to checking for 
the presence of DG. An extinctio n map covering the whole sky 
was produced by Dobash^ (1201 ll) using the Two Micron All Sky 
Survey Point Source Catalog (2MASS PSC). This map has since 
been significant ly improved to corre ct for the background-star 
intrinsic colors (Dobashi et al.L l20T2b . 

This study is carried out fo llowing the methodology adopted 
bv lPlanck collabor ationl (1201 lal) . but using extinction instead of 
FIR optical depth. The use of dust optical depth requires us to de- 
termine dust temperature and therefore make assumptions about 
the dust emissivity shape and/or mixing effects along the line-of- 
sight (LOS). On the other hand, in principle, extinction directly 
measures dust column density. In practice however, large-scale 
extinction maps derived from star counts have limited accuracy 
and suffer from a bias inherent to the stars being intermixed with 
the gas, and starlight does not sample th e whole LOS. In this 
analysi s, we use the same gas tracers as in iPlanck collaboration! 
d201 lab , so that the results can be compared directly. 

Here, we compare the spatial distribution of the DG seen 
in absorption and emission, as well as the DG masses derived 
from the two approaches. We also wish to investigate whether 
the transition between the HI and H2 phases is consistent with 
that found using FIR emission. In addition, we are interested in 
comparing our resu lts with theoretical models such as that of 
IWoifireetal.l(l2010h . 

In Section |2 we present the dataset used in this analysis, in 
Section [37X1 we then describe the method we applied to perform 
the correlations between extinction data and gas tracers. We dis- 
cuss the results in Section [4] Finally, a summary of our findings 
is provided in Section[5] 

2. Observations 

2.1. Extinction data 

We used color excess maps produced b y iDobashil (|201 ll) based 
on the 2MASS PSC dSkrutskie et all I2006I) . The maps were 
derived using a technique named "the X percentile method" 
dDobashi et all I20081 l2009t) . which is an extension of t he well- 
known NICE method introduced by lLada et all (Il994l) . While 
the standard NICE method uses the mean color of stars found 
in a cell set on the sky to measure the color excess by dust, 
the X percentile method utilizes the X percentile reddest star 
(X = 100% is the reddest) and is unaffected by the contami- 
nation of unreddened foreground stars. 

The combination of the maps of the Galactic plane and both 
north and south high Galactic latitudes cover the whole sky with 
a 1 ' grid and a varying angular resolution (from 1 ' to 12') be- 
cause of the "adaptive grid" mapping technique used (Cambresv, 
Il999l) . This technique adjusts the resolution to ensure a constant 
number of stars in the cells while measuring the color excess, 
and maintain a flat noise level all over the sky. 



Among the color excess maps presented bv IDobashil d201 ll) . 
we chose to compare the E(J - H) map measured in the range 
50 < X < 95 % (based on his equation 4) with the gas data, 
because this map is relatively more sensitive than other maps 
and has fewer defects in the regions studied in this paper (\b\ > 
10°). However, we note that the results obtained in the following 
do not change significantly when using other maps computed 
at different X values or from the other color excess available, 
E(H - K s ). 

Here, we should note that there is a slig ht but sy s temat ic off- 
set in the original color excess maps of IDobashil Q201 1 ); this 
arises from the rather ambiguous determination of the back- 
ground star colors, i.e., the mean star colors unreddened by dust, 
which is needed to determine the zero point of the color excess. 
In general, when deriving a color excess map on large scales us- 
ing the NICE method, it is virtually impossible to determine the 
background star colors precisely from the observed stars them- 
selves, because the mean intrinsic star colors should vary on the 
sky and there i s no region with out dust across the Galactic plane. 
For simplicity, IDobashil d201 ll) assumed a constant mean intrin- 
sic star color in all directions, as for ot her all-sky color excess 
maps derived in a similar way, e.g., by iRowles and Froebrich I 
(2009); this led to a systematic offset in the final color excess 
map. 

This varying offset is generally not large, but can cause 
significant problems when attempting to detect DG. To reduce 
these o ffsets, we use the new color-excess maps of Dobashi et all 
(120121) . who provided a correction to the original color-excess 
maps of IDobashil d201 ll) by determining the mean intrinsic star 
colors as a funct ion of Galac t ic coo rdinates using the Besancon 
Galactic model dRobin et all 120031) . which is one of the latest 
Galactic stellar population synthesis model. On the basis of this 
model, they generated a star catalog equivalent to the 2 MASS 
PSC, but free from any interstellar dust, and calculated how the 
mean intrinsic star colors vary across the sky when applying the 
X percentile method to the simulated star catalog. They regarded 
the resulting star color maps as the background for the color- 
excess maps derived from the 2MASS PSC. 

In addition, Dobash i et all d2012l) estimated the fraction of 
total dust along the LOS that escapes detection by the X per- 
centile method, by simulating the effect of a diffuse dust disk 
on the simulated star catalog. The method is known to under- 
estimate the diffuse dust component extending over a large re- 
gion along the LOS, but it can retrieve extin ction from individua l 
well-defined dark clouds more precisely (Doba shi et all |2009). 
As a result, they found that more than 80 % of the total dust 
along the LOS should be detected in the Galactic latitude range 
|b| > 5°, though the effect is more severe in the lower Galactic 
latitude range (e.g., by ~ 50 % at |b| =s 0°). They derived a map 
of the detection rate of the dust disk, which can be used to correct 
the color excess maps for this underestimation on large scales. 

The color excess maps used in this work were corrected for 
the background star colors as well as for the underestimation 
of the diffuse and extended dust component. However, the cor- 
rected maps might still be erroneous close to |b| 0°, especially 
around the Galactic center, because of the inhomogeneous de- 
tection limits in the 2MASS PSC and/or unknown stellar popu- 
lations not taken into account in the Besancon Model. Therefore, 
when we compare the maps with the gas data in the follow- 
ing, we restrict ourselves to |b| > 10°, where these corrections 
for the background and the underestimation are minimized. This 
also co rresponds to the region studied in IPlanck collaboration! 
d2011al) . 
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Fig. 1. Extinction map (Ay) in mag. Regions in light blue are not covered by the combination of the Dame et al ] (1200 lh and 
NANTEN CO surveys, and are not used in our analyses. The red boxes correspond to individual regions defined in Section l3.1.3l 



Finally, we converted the E(J - H) map to Ay using 

A v = 10.9£(7 - H), (1) 

w here the coeffic i ent co rresponds to the empirical reddening law 
of ICardelli etall (1 19891) for R v = 3.1. The resulting map (see 
Figure 1) has an almost constant noise level of AAy — 0.5 mag. 

2.2. Gas tracers 

T he gas tracers used in t his ana lysis are the same as described 
in iPlanck collaboration! (1201 lal) . Here we provide a summary 
of the atomic (21 cm HI emission) and molecular ( 12 CO(J=l- 
0) line) surveys. We used the LAB (Leiden/Argentine/Bonn) 
survey to trace the atomic gas. T his survey is a comb i nation 
of the Leiden/Dwingeloo survey bv lHartmann & Burtonl (1 19971) 
with sky observations above -30° of Galactic latitude (at a 
36' angular resolution), and the IAR (Institute) Argentino de 
Radioastronomia) survey (lArnal et al.U2000llBaiaia et al.l 12005) 
of the Southern sky at latitudes below -25° (at a 30' angular res- 
olution). The LAB data were integrated in the velocity range 
-400< Vlsr<400 km s _1 . Assuming that the gas is optically thin, 
we deduced the hydrogen column density from the integrated 
intensity of the HI emission (Whi) using 



Nm = XhiWhi, 



where Xhi is the HI integrated intensity to column-density con- 
version factor. Thi s factor is tak en to be equal to 1.82xl0 18 
H/cm 2 /(Kkms" 1 ) dSpitzerifl978h . 

For the molecular gas, we used the combination of three 
12 CO(J=1-0) line surveys: 

- The iDame et"al .] (1200 lh survey, in the Galactic plane, ob- 
tained with both the CfA telescope in the north (at an angular 
resolution of 8.4') and CfA-Chile telescope in the south (at 
an angular resolution of 8.8'). The integrated intensity map 
was derived by integrating the vel ocity range w here the CO 
emission is significantly detected (Dama. 1201 lh . 



- The unpublished high latitude survey obtain ed with the 
CfA t elescope, still observing the northern sky (Dam e et"aQ 
|2012|) . The data cube was integrated over 10-20 velocity 
channels. 

- The NANTEN survey obtained from Chil e for the intermedi- 
ate Galactic latitudes not covered by the IDame et aL I (120011) 
survey, at a 2.6' angular resolution. This survey is still un- 
published, but a full descripti on of the NANTEN telescope 
can be found for instance in lFukui et al.l dl999h . The total in- 
tensity map was obtained by integrating the data cube over 
the whole velocity range. 

Each survey was smoothed to a common resolution of 8.8' 
using the convolution by a Gaussian kernel. As quoted in 
Planck collaboration (201 la), the NANTEN data appear to have 
24% higher values than those obtained with the CfA tele- 
scop es. This discrepancy is still poorly understood and follow- 
ing |Planckcon*oration| (12011 ah . the NANTEN data were de- 
creased by 24%, before mer ging with the CfA survey, to bring 
them to the absolute scale of Da me et alj d2001l). The total map 
used in this analysis is shown in Planck collaboration! (1201 1 ah . 
Figure 1 . The CO integrated intensity (Wco) to column density 
conversion factor is derived from the relation 



(2) N H2 = XcoWco- 



(3) 



The value of the CO conversion factor (Xco) is still debated. 
This factor is also expected to vary over the sky. For this reason, 
we made no assumption and derived Xco instead from correla- 
tions (see Section lXTT i. 

The contribution of emission in the ionized medium is 
neglected in this study. This assumption is reasonable since we 
do not consider the Galactic plane emission, which includes 
most of the ionized medium emission caused by HII regions. 
However, its possible small contribution at high latitudes, which 
is correlated neither with atomic nor molecular-traced gas, can 
be accounted for in the constant term (A y ) when performing the 
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correlations (see Equation|4]l. 

This work is done at the resolution of the HI data, i. e. 
36'. All data are projected using the HEALPix pixelization 
scheme (Hierarchical Equal Area isoLatitude Pixelization£| with 
nside=256, corresponding to a pixel size of 13.7'. The descrip- 
tion of the projection method used is given in AppendixlAl 



3. Extinction/gas correlation 

3.1. (A v /N H ) ref andX C o determination 

In lPlanck collaboration! (|201 lal) . the reference value of the ratio 

of the dust optical depth to the gas column densities ((j^j-) ) 
was determined from the correlation between t and the gas trac- 
ers at low gas-column density. We followed the same type of 



analysis, and modeled the extinction (A™ ) in the portion of the 
sky covered by the infrared, atomic, and CO-traced surveys, for 
|b|>10°, as 



i mod 



+ 2X co W co 



(4) 



(A \ re f 
^ 1 is the reference value of the ratio of visible ex- 
tinction to gas column density, and Ay is a constant, which can 
account for ionized gas and/or the offset in the extinction map. 
However, the Ay map is particularly noisy at low col- 

(A \ re ^ 
j^j diffi- 
cult. We therefore a d opted a slightly different approach from 
iPlanck col laboration! d201 lab . Our analysis was done in two 
steps. We first selected regions wi th a limited amount o f DG by 
identifying regions in Figure 8 of Planck collaboration! (1201 lal) 
with DG column densities N* such that N^/Nh < 0.05 and we 
also imposed a low CO content (Wco < 0.2 K km s" 1 ). For this 
region, which we refer to as "no DG", we derived the best-fit 

parameters (j^j and A° using Eq.|4l but did not attempt to 
derive Xco in this first step because the considered regions have 
little CO emission. Second, we repeated the same analysis in- 
cluding all pixels of the considered region included in the gas 
surveys (we refered to this region as "all"), while imposing the 

(jifa) and A y values derived in the first step and determined the 
best-fit value for Xco- We also computed the mass of DG from 
the difference between the best correlation and the data at that 
stage. The procedure above was applied to both the whole high 
Galactic latitude sky (see Sec. 13.1.11 ) and a set of individual re- 
gions (see Sec. 13.1.31 . and in all cases was restricted to |b| >10°. 

Results of the correlations are presented in Table [T] The best- 
fit parameters were obtained from x 1 minimization. Parameter 
uncertainties were derived from the difference between the min- 
imum and maximum values of the parameters contained in in- 
terval Ay 2 , corresponding to a confidence level of 68%. We note 
that, for the purpose of illustrating the Ay-gas correlations (see 
Fig. |2] and |3), we also carried the correlation in a region with 
N H /N H > 0.05, which we refered to as the "DG" region. 

3.1.1. High Galactic latitude regions (|b| > 10°) 

The correlation was first established across the whole high 
Galactic latitude sky (|b| > 10°) characterizing the solar neigh- 
borhood. Results are displayed in Figure |2] (top panels). We 



/ a \ref 

note that the linear fits (slope=(^j ) shown in the figures 
(red lines) were constrained using the whole data in the "no 
DG" region, and not the binned version of the data shown by 
yellow dots. The corresponding lines are curved as plotted on 
a log-log scale, because of the constant term. The investiga- 
tion of this "no DG" region leads to an (Av/Nh) ratio of 
6.53xl0" 22 mag. cm 2 . Thi s value is close to t he reference value 
of 5.34x 10" 22 obtained bv lBohlin et all d!978l) . which is relevant 
to the solar neighborhood and the Galactic plane. We found that 
the optimal Xco is close to 1.67 X 1 20 H 2 cm' 2 /(K km s~ 1 ) . Thi s 
value is lower than the findings of Planck collaboration! (1201 lal) . 
who derived an averaged Xco=2.54xl0 20 H2cm" 2 /(Kkms _1 ), 
but in better agreem ent with the Galac t ic av erage of 1.9 xlO 20 
H 2 cm _2 /(Kkms _1 ) dStrong & Mattoxl [19961) and the valu e of 
1 . 8 x 1 20 H 2 cm 2 /(K km s" 1 ) derived bv lDame et ail d200ll) for 
|b|>5°. 

Looking at Figure [2] one can discern a different behavior in 
the Ay versus Nh plots (with Nh = N™+2XcoWcoX when com- 
paring the "DG" and "no DG" regions at intermediate and high 
Ay. In the first case, looking at the binned data points the cor- 
relation shows an excess between Nh - 4 X 10 20 crrT 2 (Ay =0.2 
mag) and Nh - 3 X 10 21 cm" 2 (Ay=1.5 mag) compared to a 
linear fit, whereas we do not in the second case. We note that, 
owing to the low signal-to-noise ratio of the extinction data at 
low column densities, the location of these transitions were de- 
termined by eye here, while they were inferred from a fit in 
IPlanck col laboration! (1201 lal) . These results imply that there is 
an additional gas phase, with a spatial distri bution that is similar 
to tha t determined using the Planck data (IPlanck collaboration, 
1201 lal) . For comparison, these authors estimated that the excess 
appears at Ay ranging from 0.4 to 2.5 m ag. However, results ob- 
tained bv lPlanck collaboration! d201 lbl) . by analyzing the diffuse 
ISM and Galactic halo, revealed some residual emission in the 
infrared/HI correlation, between N H - 3 x 10 20 cm" 2 (A v ^0.15 
mag) and Nh - 4 x 10 21 cm' 2 (Ay ^2.0 mag). Our results are 
in-between these two findings based on Planck data. 

The column density of the dark component (N*) can be com- 
puted as the difference between the observed and the best-fit 
(A™ od ) extinction per unit column density, over the entire region 



4„ _ A" wd 
_ Ay Ay 

H (A v /N H y<f 



(5) 



The total mass of hydrogen within the DG component (M^), the 
atomic (M™), and C O-traced molecular ga s (M^ °), is estimated 
in the same way as in IPlanck collaboration! (1201 lal) 



(6) 



http://lambda.gsfc.nasa.gov/ 



where D is the distance to the gas, mn the hydrogen atom 
mass, a nd Q p ; x the pixel solid angle. As in IPlanck collaboration! 
d201 lal) . we assumed the same distance to all gas components. 
We obtained M^/M™=0. 19+0.25 and M|/Mg°=1.64±2.20. 
These values are expected to be both overestimated and underes- 
timated, respectively, since the scale-height of the HI component 
is larger than that of the molecular one, with a probable inter- 
mediate scale-height for the DG component. The formal uncer- 
tainty in the dark gas mass is large and essentially results from 
the assumed uncertainty in AAy /Ay and from our assumption 
that this uncertainty is an absolute error and therefore does not 
scale down with the number of independent measurements. This 
leads to a large uncertainty in the DG mass when most pixels in 
a region have low Ay values, but it is likely that our assumption 
about the nature of the errors overestimates the uncertainty. 
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Outer Galaxy 


5.46+0.03 


4.68±0.75 


1.67+0.10 
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0.43 


0.3* 
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3.41+0.11 


7.47±0.29 


2.84±0.92 


0.45±0.12 


0.66±0.18 


0.40 


0.1* 


Aquila-Ophiuchus 


6.99+0.18 


11.78+0.53 


3.46±0.91 


0.59+0.06 


1.65+0.18 
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0.6* 



Table 1. Derived parameters and their l-cr uncertainties, computed over different regions, for |b| > 10°. 



" From Planck collaboration d2011ah 
* From Grenier et al. (2005) 
c From Ubdo. et al.1 d201(Hl 




Fig. 2. Correlation plots between extinction data (Ay) and the total gas-column density (Nh — N™ + 2XcoWco) excluding the 
Galactic plane (|b| > 10°): high latitude sky (top panels), inner Galaxy (middle panels) and outer Galaxy (b ottom panels). Left, 
middle, and right panels correspond respectively, to regions where no DG is detected ("no DG") according to iPlanck collaboration 
j20TTah . where strong DG is detected ("DG"), and the entire region ("All") (see Sec. 13. U . The color scale represents the density 
of pixels on a log scale. The yellow dots show the correlation binned in Nh- The red lines show the best-fit linear correlation, 
constrained for the whole data in the "no DG" region. The short-dashed and long-dashed lines indicate the position of a departure 
from a linear fit, and the position where the correlation recovers a linear behavior, respectively. 



Our DG mass is lower than found in Planck collabo ration! 
( 1201 lab , who deduced M^/M™ 0.28 ± 0.03 and M£/M£° 
1.18 ± 0.01, but consistent with it, taking our large uncertainties 
into account. With the use of their Xco value, we would obtain 

M h/ M h° ~ °- 73 - ° ur results lead t0 a DG mass e q ual t0 17% 
of the total observed mass, close to the value of 22% derived 



bv lPlanck collaboration! (1201 lab . Our study using extinction data 
therefore reveals slightly less DG mass than the same analysis 
based on FIR data. The difference between the two analyses is 
however insignificant owing to large uncertainties for this region 
of the sky. 
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Fig. 3. Correlation plots between extinction data (Ay) and the total gas column density (Nh = N™+2XcoWco) in various regions ex- 
cluding the Galactic plane (|b| > 10°), from top to bottom the Cepheus-Polaris, Taurus, Orion, and Aquila-Ophi uchus regions. Left, 
middle, and right panels correspond respectively, to regions where no DG is detected ("no DG") according to iPlanck col laboration 
la) , where strong DG is detected ("DG"), and the entire region ("All") (see Sec. 13. U . The color scale represents the density of 
pixels on a log scale. The description of the yellow dots, red, and dashed lines is given in the caption of Figure|2] 



3.1.2. Inner/outer Galaxy 

We carried out the same analysis for the inner (|1| < 70°) and 
outer (|1| > 70°) Galaxy, at high Galactic latitudes (|b| > 10°). 
A difference appears in the (Ay/NH) ref ratio, with a higher 

value being found in the inner Galaxy ((^Y = 8.72 x 1(T 22 

mag.crrT 2 ) than the outer Galaxy ((^)' ef = 5.46 x 1(T 22 
mag.cirT 2 ). This result could illustrate the true large-scale varia- 
tions in the dust abundance and reflect the large-scale metallicity 
gradient in our Galaxy. 

From Figure [2] one can see that the excess is more promi- 
nent in the inner Galaxy than the outer Galaxy. In addition, we 
note that the outer Galaxy appears to be more noisy at low gas 
column densities. This is likely because the outer Galaxy hosts 
fewer stars than the inner Galaxy, inducing larger uncertainties 



in the resulting extinction map. The reconciliation between the 
data and the linear fit, corresponding to the H2-to-CO transition 
in the case of molecular DG, seems to occur at Nh — 3-3.5xl0 21 
cm -2 (Ay —1.5-1.75 mag) in both the inner and the outer Galaxy. 



In term of the DG mass, each part of the Galaxy has a dif- 
ferent behavior. The inner regions contain a more enhanced DG 
component, with ratios M^/M™ and M H /M H ° of about 0.24 and 
2.43, respectively. In contrast, the DG mass in the outer Galaxy 
represents only 1 1% of the gas in the atomic phase and 74% of 
that in the CO-traced phase. This result indicates that the addi- 
tional gas phase is essentially concentrated in the inner Galaxy. 
In a similar way, the Xco factor appears to be larger in the inner 
Galaxy (X co = 2.28 x 10 20 H 2 cm^AKkms- 1 )) than its outer 
parts (X co = 1.67 x 10 20 Hzcm^/^kms- 1 )). 
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3.1.3. Individual regions 

Performing the analysis on individual regions requires the avail- 
ability of enough pixels with and without DG, taking into ac- 
count the limited extent of the CO survey. We were able to ex- 
plore the Ay-Nn correlation in the following four regions: 

- Taurus (200° < I < 160°; -35° < b < -10°); 

- Orion (240° < I < 200°; -30° < b < -10°); 

- Cepheus-Polaris (135° < I < 90°; 10° < b < 40°); 

- Aquila-Ophiuchus (|Z| < 20°; 10° < b < 48°); 

The results of the correlations for individual regions are pre- 
sented in Figure [3] Each environment has different dust proper- 
ties, with i^f ranging between 3.41xlCT 22 and 6.99 xl(T 22 

mag cm 2 . Moreover, these regions are of great interest because 
they have distinctive gas properties. The Xco factor in Orion and 
Aquila-Ophiuchus are larger than in the other regions. In addi- 
tion, the former regions have a substantially large DG mass ratio 
relative to the atomic hydrogen, with ratios of 0.45 and 0.59, 
respectively. These large amounts of DG in Orion and Aquila- 
Ophiuchus are visible in Figure [3] as a clear, large excess in the 
Ay-Ne plot. We checked that the excess was not due to con- 
tamination by young stellar objects (YSOs). The color excess 
maps may indeed overestimate the true Ay if there is significant 
contamination by young stellar objects (YSOs). In these two re- 
gions, we therefore also performed the correlation between ex- 
tinction maps in the J band (Aj) and the gas tracers, because 
YSOs would induce an underestimation in the extinction maps 
(as opposed to the color excess maps). These results indicate that 
the large amount of DG is not caused by the presence of YSOs. 
iGrenieret ai1 (l2005) obtained a ratio M^/M™ that is even larger 
than our results, in the Aquila-Ophiuchus-Libra region, with a 
ratio of 1 . In contrast, they derived a value of 0. 14 in Orion. They 
found a ratio of 0.33 for the Taurus-Perseus-Triangulum region, 
close to our value of 0.42 in the Taurus region. We no t e how ever 
that the individual regions defined in Gren ier et al. 1 (12001 are 
not exactly the same as the ones used in this paper. The Ay-Nn 
plot for Cepheus-Polaris does not display any noticeable excess. 
When comparing the DG mass with the CO one in the different 
clouds, we observe significant variations. We note values rang- 
ing from 0.17 in the Cepheus-Polaris region to 1.65 in Aquila- 
Ophiuchus, with values of 0.66 and 0.76 in Orion and Taurus. 
These ratios are however highly dependent on the value of the 
Xco factor. The low uncertainty in the dark mass of the Taurus 
region is caused by the relatively large Ay values. The absolute 
uncertainty AAy/Ay is larger in Cepheus-Polaris which corre- 
sponds to a large number of pixels at low Ay, inducing a consid- 
erable mass uncertainty. 

The H2-to-CO transition, which we determined by eye, ap- 
pears at Nh — 3.5 x 10 21 cirT 2 (Ay =* 1.75 mag) in Taurus, 
N H =* 6 x 10 21 cm" 2 (A v ^3.0 mag) in Orion, Nh — 3.5 x 10 21 
crrT 2 (A v =a 1.75 mag) in Cepheus-Polaris, and N H - 4 x 10 21 
cirT 2 (Ay =i 2.0 mag) in Aquila-Ophiuchus. This transition 
seems to vary with environment. The position of the HI-to-H2 
transition is uncertain because of the noise in the extinction data 
at low gas-column densities. For this reason, we do not discuss 
its value further in this paper. 

4. Discussion 

The two quantities r and Ay are in principle proportional. 
However, the estimate of the dust optical depth t requires knowl- 
edge of the dust equilibrium temperature, requiring an additional 



step in the study of the dust/gas ratio, and additional assumptions 
that could bias the results. For instance, a single dust tempera- 
ture is often assumed along the LOS. In addition, the temperature 
estimate depends on the assumption made about the emissivity 
spectral index f3. In that sense, the use of extinction data is more 
straightforward. However, extinction data suffer from several de- 
fects, especially when producing the color excess map. Even if 
the X percentile method used to generate the map is a promising 
method, the resulting map can still be biased by missing extinc- 
tion, inducing variations in the detection rate. The detection rate, 
i.e., the fraction of the detected Ay to the true total Ay integrated 
along the LOS, could vary as a function of both the Galactic co- 
ordinates and the distance to the clouds. To determine how much 
extinction we should miss in the adopted Ay map, we have car- 
ried out a test by applying the X percentile method to some ar- 
tificial clouds set in the distribution o f stars generated using the 
Besan^on Model (Robin et aUl2003l) . We first prepared a spher- 
ical cloud having a Gaussian density distribution with a size of 6 
pc (at FWHM) and a total extinction along the LOS of A v = 10 
mag at the center of the cloud. We then located the cloud at re- 
alistic coordinates and distances of some nearby clouds in the 
simulated star distribution , e.g., at (l,b)=(174.0°, -13.5°) and at 
D=150 pc for the Taurus cloud, and performed the X percentile 
method exactly in the same way as iDobashU d201 lb did using 
the 2MASS PSC. Our results indicate that most of the individ- 
ual clouds analyzed here are well-detected with a detection rate 
larger than 90% (i.e., the model cloud with Ay = 10 mag is de- 
tected as Ay > 9 mag), except for Orion whose detection rate 
is ^76%. General results show that nearby clouds (D<500 pc), 
which represent most of the clouds at |b| > 10° in the Ay map, 
are clearly detected with a detection rate > 80% for all direc- 
tions. The low (Ay/NH) ref ratio derived in Orion could be the 
result of underestimating Ay in this region. Apart from this re- 
gion, the (Ay /Nn) lef ratios i n Table [Jappear eq ual or larger than 
the reference value given bv lBohlin et al.l(ll978l) . confirming that 
the X percentile method does not miss a substantial amount of 
extinction. 

The determination of the Xco conversion factor dif- 
fers from one analysis to another. Our average value of 
1.67xl0 20 H2 cm- 2 /(Kk ms -1 ) is close to the Galac tic average 
(IStrong & MattoxL Il996) and the value derived by iDame et al.l 
(2001) for |b| > 5°. A decrease in Xco is observed from the in- 
ner to the outer high Galactic latitude sky. Most of the clouds 
observed in the Galactic extinction map at |b| > 10° are nearby 
local clouds, at a distance of about 200-500 pc, so we do not at- 
tribut e the variations to the metallicity gr adient ([R olleston et all 
2000), but instead to local variations. Arimo to et al.l (11996). 
using CO data and virial masses, derived an average value 
near Sun of 2.8xl0 20 H2 crrT 2 /(Kkm s _1 ), with variations rang- 
ing from 2.09xl0 20 to 3.74xl0 20 H? cm" 2 / ( K fans " 1 ) with 
increasing Galactocentric radius. lAbdo. etafl (1201 0b discov- 
ered a similar behavior, but with significant variations rang- 
ing from 0.87xl0 20 H 2 cm- 2 /(Kkms 1 ) in the Gould Belt to 
1.9xl0 20 H2 cm~ 2 /(Kkm s" 1 ) in the Perseus arm. Finally, the 
iPlanck col laboration (2 01 lal) estimated an average value of 
2.54X10 20 H2cnr 2 /(Kkms _1 ) in the solar neighborhood. We 
also note that some of the studies included the dark component 
in their derivation of Xco- 

In their model to estimate the DG mass, WolfireetaL|([2010) 
predicted an HI-to-H2 transition located at Ay 0.2 mag, which 
agrees with our findings. This result is however slightly lower 
than the value derived using FIR data (Ay 0.4 mag) following 
Planck collaboration (201 la). However, in the latter case the HI- 
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to-H2 transition was determined by y 2 minimization, whereas in 
our study it was done by eye, which can induce large uncertain- 
ties. 

Moreover, the model predictions suggest a constant fraction 
of the molecular mass in the dark component of about 0.3 for 
an average extinction Ay around 8 mag. In cases of decreasing 
extinction, they found an increase in the dark mass fraction. In 
the framework of their model, this would indicate that a larger 
fraction of molecular gas is located outside the CO region. This 
fraction (fDG) is computed as 

Mf 

fDG - 7TF7J- 

Values of foG for each region are provided in Table [1] Over the 
entire high Galactic latitude sky, we derived foG = 0.62, with a 
larger value in the inner Galaxy (foG = 0.71) with respect to the 
outer Galaxy (foG = 0-43 ). For comparison, using the s ame defi- 
nition of this quantity, the lPlanckc ollaborationl d201 la ) obtained 
a value of fpG —0.55 for the solar neighborhood. The Taurus 
and Aquila-Ophiuchus regions have a DG fraction of 0.43 and 
0.62, respectively. These values are in g ood agreement with the 
results obtained bv lGrenier et al.l (120051) (see Table [1} in equiv- 
alent regions, with fractions of 0.3 and 0.6, for the Taurus and 
Aquila-Ophiuchus-Libra regions. These authors derived fuc as 
low as 0. 1 for Cep heus-Cassiopeia-Polaris and Orion, whereas 
lAbdo. et aDJ20 10) deduced a value of 0.30 for Cepheus. We ob- 
tained foG equal to 0.15 and 0.40 in Cepheus-Polaris and Orion. 

All these studies illustrate the existence of variations in the 
dark-gas mass fraction. In our study, we have excluded the 
Galactic plane, where the most massive and bright molecular 
clouds are located. At the 36' angular resolution, most of the 
pixels considered in our analysis have extinctions lower than 10 
mag. I n this case, fpr, va lues larger than 0.3 are expected fol- 
lowing and this is indeed what is observed 
here. 

Dust in the form of aggregates is likely to be more emissive 
than isolated grains, but aggregation is not expect ed to substan- 
tially affect absorption properties in the visible (iKohler et al.L 
1201 ll) . The detection of significant departure from linearity be- 
tween extinction data and gas column densities, especially in 
Orion or Aquila-Ophiuchus, which exhibit the largest excess 
in Figure [3] indicates that dust emissivity changes induced by 
grain aggregation is an unlikely explanation of the observed ex- 
cess. As a consequence, grain coagulation can only be respon- 
sible for a small amount of the observed excess, and may ex- 
plain some of the differences in the studies based on FIR data 
dPlanck collaborationi 1201 lal) and visible extinction data (this 
w ork). Moreover, our results agree with the model prediction 
of IWolfire et al.l (1201 Ol) . favoring the hypothesis of pure H2 gas 
(without CO molecules), surrounding the CO regions. 

5. Summary 

Using a recently revised all-sky extinction map, we have exam- 
ined the correlation between the extinction and gas column den- 
sities derived from HI and CO observations. Our results can be 
summarized as follows: 

- We have measured an excess of extinction at intermedi- 
ate column densities, relative to a linear-fit between Ay 
and observable gas. This excess is observed over the whole 
high Galactic latitude sky, with a predominant amount being 
found toward the inner regions (|1| < 70°), and in particular 



in Aquila-Ophiuchus. This result confirms the recent detec- 
tion of DG in the Planck data and implies that the effect in 
the FIR is not due to changes in the dust optical properties 
caused by dust aggregation. 

- We have derived an average dust extinction to gas ratio of 
(A v /N H ) ref = 6.53 x 10~ 22 mag cirT 2 in the solar neighbor- 
hood, for an average Xco = 1.67 x 10 20 H2crrT 2 /(K km s _1 ). 
The inner/outer Galaxy exhibits a higher/lower (Ay /NH) ref 
ratio (8.72xl0"~ 22 /5.46xl0~ 22 mag cm -2 ), associated with a 
higher/lower X co factor (2.28xl0 20 /1.67xl0 20 H 2 cm" 2 /(K 
km s _1 ). A significantly larger value, more than twice the 
value derived in the solar neighborhood is observed in the 
Aquila-Ophiuchus region. In contrast, the lowest Xco value 
is found in Cepheus-Polaris. 

- The results of our analysis agree with the theoretical predic- 
tions of modeling of the dark H2 gas obtained to help explain 
the observed departure from linear-fits in the correlations . 

- The DG mass derived is 19% of the atomic mass and 164% 
of that in the CO-traced gas in the solar neighborhood. It can 
reach up to 59% of the atomic gas in the Aquila-Ophiuchus 
region. Our DG mass estimates are slightly lower than that 
derived from FIR Planck data, but the difference is insignifi- 
cant. 

- We have estimated the fraction of the molecular mass in 
the dark component (fDG) and found an average value of 
0.62 in the solar neighborhood, with variations going from 
0.43 to 0.71 in the outer and inner Galaxy. We have de- 
rived f D o=0.15, 0.43, 0.40, and 0.62 in the Cepheus-Polaris, 
Taurus, Orion, and Aquila-Ophiuchus regions, respectively. 

- The HI-to-H2 and H2-to-CO transitions appear for Ay =* 0.2 
and Ay =s 1 .5 mag in the solar neighborhood, with variations 
from regions to regions. 
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Appendix A: WCS to HEALPix ancillary data 
transformation 

The data used in this paper were transformed from their na- 
tive WCS (world coordinate system) local projection into the 
HEALPix (Hierarchical Equal Area isoLatitude Pixelization) 
all-sky pixelization using the method described here in this 
Appendix. This method has also been used to produce 
HEALPix maps of a larger set of ancillary data, particularly in 
the context of the analysis of the Planck data0. We note that the 
CO map used in this paper is not yet freely available, since the 
corresponding data have not been made public. 

The HEALPix format allows us to store ancillary data on a 
single grid scheme on the sphere. This is advantageous since dif- 
ferent data can then be compared on a pixel by pixel basis, with- 
out the need to project the data to a common grid. Ancillary data 
are available in the WCS convention and are usually stored as 
Flexible Image Transport System (FITS ) files. In the following, 
we refer to the ancillary data as FITS maps. It is usually imprac- 
tical to return to the raw ancillary data and directly reprocess the 
maps into the HEALPix pixelization. It is therefore necessary to 
project the FITS format data onto the HEALPix pixels. As the 
FITS and HEALPix pixels do not match each other in position, 
size, and shape, this requires some kind of interpolation, which 
must be done with the minimal loss of information and without 
altering the photometry of the original ancillary data. Here, we 
use a mosaicking method where we compute explicitly the sur- 
face of the pixel intersections, and use these values as weights to 
construct the HEALPix map. 

The HEALPix pixels projected onto a local FITS map are 
shown in Figure I A. 1 1 The calculations are performed for a 
HEALPix pixel size (given by the N s jd e HEALPix parameter) so 
that the corresponding HEALPix pixel size matches the Shannon 
criterion for the angular resolution of the data considered. This 
ensures that no spatial information is lost in the conversion 
and that the corresponding HEALPix and the FITS pixels have 
similar sizes. For a given FITS projection (i.e. from the as- 
trometry information contained in the FITS header), we first 
identify the HEALPix pixels that intercept the sky area cov- 
ered by the FITS map. For each of those HEALPix pixels, we 
identify the FITS pixels with a non-zero intersection. We then 
compute the surface fraction of the FITS pixels intersecting the 
HEALPix pixel Sjjj, where "if" denotes the FITS pixel number 
and "ih" denotes the HEALPix pixel number. 

The calculation is performed in the pixel space of the 
FITS map shown in Figure I A. 1 1 We first compute the sky co- 
ordinates of the four corners of the considered HEALPix pixel. 
We then transform these coordinates into 2D pixel numbers (i,j) 
of the FITS image using standard routines and the astrometry 
information contained in the FITS header of the considered im- 
age. We then assume that the frontier of the HEALPix pixel is a 

2 The ancillary data can be accessed on 
http : //www . cesr . fr/-bernard/ Ancillary/ 
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Fig. A.l. Geometry of the HEALPix and FITS pixels. 

straight line in the FITS image pixel frame. Although this is not 
exactly true, this is a very small approximation when the size of 
the HEALPix pixel is smaller or on the order of the FITS pixel 
size, which is always the case here. We then compute the sur- 
face of the corresponding intersection polygon, which has been 
normalized to that of the HEALPix pixel surface SjjjB 

The intersection geometry informations (if and S]^) are 
stored in a table containing one entry per HEALPix pixel 
that has non-zero intersection with the considered FITS map. 
Calculations are carried and stored separately for each 
FITS map. We note that the above table can be easily inverted to 
a table with one entry per FITS pixel, giving the HEALPix pixel 
numbers intercepting it (ih) and the corresponding surface frac- 
tion Sjj?, which in turn can be used to perform reverse calcula- 
tions projecting HEALPix maps onto local WCS projection. 

The above intersection fractions are then used to compute 
HEALPix ancillary data value djh, given the FITS data dif as 

d ih = J] S "f xdi f> (A -V 

if 

where the summation is carried over all FITS pixels intersecting 
a HEALPix pixel ih. 

When the data is composed of a collection of individual 
FITS maps, the calculation are performed for each map sepa- 
rately, while maintaining a map of the total weight, which al- 
lows us to evaluate average values in sky regions where indi- 
vidual maps overlap. We note that virtually all WCS projection 
types can be used, including the sixcube projection used for the 
COBE and FIRAS data. 



3 This calculation is performed using the polyfillaa routine (see 
http : //tir . astro . utoledo . edu/ j dsmi th/code/idl . php ) opti- 
mized for fast clipping of polygons against a pixel grid. 



